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We report the first measurements of the kurtosis (k), skewness (S) and variance (a 2 ) of net- 
proton multiplicity (iV p — Np) distributions at midrapidity for Au+Au collisions at ^/snn = 19.6, 
62.4, and 200 GcV corresponding to baryon chemical potentials (/j.b) between 200 - 20 MeV. Our 
measurements of the products Ka 2 and Sa, which can be related to theoretical calculations sensitive 
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to baryon number susceptibilities and long range correlations, are constant as functions of collision 
centrality. We compare these products with results from lattice QCD and various models without 
a critical point and study the -^/snn dependence of na 2 . From the measurements at the three beam 
energies, we find no evidence for a critical point in the QCD phase diagram for (ib below 200 MeV. 

PACS numbers: 25.75.Gz,12.38.Mh,21.65.Qr,25.75.-q,25.75.Nq 



One of the major goals of the heavy-ion collision pro- 
gram is to explore the QCD phase diagram [1]. Finite 
temperature lattice QCD calculations Q at baryon chem- 
ical potential /ib = suggest a cross-over above a critical 
temperature (T c ) ~ 170 - 190 MeV Q from a system 
with hadronic degrees of freedom to a system where the 
relevant degrees of freedom are quarks and gluons. Sev- 
eral QCD based calculations (see e.g find the quark- 
hadron phase transition to be first order at large p,B ■ The 
point in the QCD phase plane (T vs. /ib) where the first 
order phase transition ends is the QCD Critical Point 
(CP) Hi]. Attempts are being made to locate the CP 
both experimentally and theoretically [7]. Current the- 
oretical calculations are highly uncertain about location 
of the CP. Lattice QCD calculations at finite hb face nu- 
merical challenges in computing. The experimental plan 
is to vary the center of mass energy (^/snn) of heavy- ion 
collisions to scan the phase plane [8J and at each energy, 
search for signatures of the CP that could survive the 
time evolution of the system Q . 

In a static, infinite medium, the correlation length (£) 
diverges at the CP. £ is related to various moments of the 
distributions of conserved quantities such as net-baryons, 
net-charge, and net-strangeness jlu| . Typically variances 
(cr 2 = ( (AN) 2 )] AN = N - M; M is the mean) of these 
distributions are related to £ as cr 2 ~ £ 2 [llj]. Finite size 
and time effects in heavy-ion collisions put constraints 
on the values of £. A theoretical calculation suggests £ 
« 2-3 fm for heavy-ion collisions 12[. It was recently 



shown that higher moments of distributions of conserved 
quantities, measuring deviations from a Gaussian, have a 
sensitivity to CP fluctuations that is better than that of 
cr 2 , due to a stronger dependence on £ [Hj]. The numer- 
ators in skewness (S = ((AN) 3 ) /a 3 ) goes as £ 4 ' 5 and 
kurtosis (k — [((AN) 4 ) /a 4 } - 3) goes as £ 7 . A crossing 
of the phase boundary can manifest itself by a change of 
sign of S as a function of energy density [lj, [3] • 

Lattice calculations and QCD-based models show 
that moments of net-baryon distributions are related to 

baryon number (A7Vb) susceptibilities (xb = ^ A y^ ] 
V is the volume) [15[ . The product no 1 , related to the 

ratio of fourth order (Xb ) *° secon d order (Xb ) sus " 
ccptibilities, shows a large deviation from unity near the 
CP |15j . Experimentally measuring event-by-event net- 
baryon numbers is difficult. However, the net-proton 
multiplicity (N p — Np = A7V P ) distribution is measur- 
able. Theoretical calculations have shown that AN P fluc- 
tuations reflect the singularity of the charge and baryon 



number susceptibility as expected at the CP [16]. Non- 
CP model calculations (discussed later in the paper) show 
that the inclusion of other baryons does not add to the 
sensitivity of the observable. This letter reports the first 
measurement of higher moments of the AN p distribu- 
tions from Au+Au collisions to search for signatures of 
the CP. 
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FIG. 1: (Color online) AN P multiplicity distribution in 
Au+Au collisions at ^/snn = 200 GeV for various collision 
centralities at midrapidity (| y |< 0.5). The statistical errors 
are shown. 



The data presented in the paper are obtained using 
the Time Projection Chamber (TPC) of the Solenoidal 
Tracker at RHIC (STAR) 0. The event- by- event pro- 
ton (N p ) and anti-proton (Np) multiplicities are mea- 
sured for Au+Au minimum bias events at ^/snn = 19.6, 
62.4, and 200 GeV for collisions occurring within 30 cm 
of the TPC center along the beam line. The numbers of 
events analyzed are 4xl0 4 , 5xl0 6 , and 8xl0 6 for ^/snn 
= 19.6, 62.4, and 200 GeV, respectively. Centrality selec- 
tion utilized the uncorrected charged particle multiplicity 
within pseudorapidity | rj \ < 0.5, measured by the TPC. 
For each centrality, the average numbers of participants 
((-Wpart)) are obtained by Glauber model calculations. 
The AN p measurements are carried out at midrapidity 
(| y | < 0.5) in the range 0.4 < p T < 0.8 GeV/c. Ioniza- 
tion energy loss (dE/dx) of charged particles in the TPC 
was used to identify the inclusive p(p) [l8j]. To suppress 
the contamination from secondary protons, we required 
each p(p) track to have a minimum pt of 0.4 GeV/c and 
a distance of closest approach (DCA) to the primary ver- 
tex of less than 1 cm [18(. The pr range used includes 
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approximately 35-40% of the total p + p multiplicity at 
midrapidity. AA p was not corrected for reconstruction 
efficiency. Typical AiV p distributions from 70-80%, 30- 
40%, and 0-5% Au+Au collision centralities are shown in 
Fig.ffl 
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FIG. 2: (Color online) Centrality dependence of moments of 
AiVp distributions for Au+Au collisions at ^/snn = 1 9.6, 62.4, 
and 200 GeV. The lines are the expected values from the 
central limit theorem. Error bars are statistical and caps are 
systematic errors. 

The four moments (AT, a, S, and k) which describe 
the shape of the AA p distributions at various collision 
energies are plotted as a function of (A part ) in Fig. [2] 
The typical statistical errors on a, S , and k for central 
Au+Au collisions at 200 GeV are 0.2%, 1 1 % and 16% re- 
spectively. The M shows a linear variation with (N part ) 
and increases as yj snn decreases, in accordance with the 
energy and centrality dependence of baryon transport Q . 
The variation of M within a centrality bin has been taken 
into account in higher moment calculations. The a in- 
creases with (iVp ar t). The values are similar for three 
beam energies studied. The S is positive and decreases 
as (-/V part ) increases for a given collision energy. The val- 
ues also decrease as ^/snn increases. This indicates that 
the distributions become symmetric for more central col- 
lisions and for higher beam energies. The k decreases as 
(iV par t) increases, but is similar for all three ^/snn stud- 
ied. 

Experimentally it is difficult to correct such observ- 
ables for the particle reconstruction efficiency on an 
event- by- event basis. Construction of observables in- 
dependent of the efficiency, such as factorial moments, 
leads to loss of one-to-one correspondence with higher 
moments (l9| . and significant difficulty in comparing to 
theoretical expectations. We have investigated the ef- 
fects of the detector and track reconstruction efficiencies 



by comparing the moments of the AA p distribution us- 
ing the events from a heavy-ion event generator model 
HIJING (ver.f.35) [2(| and the moments of the recon- 
structed AA p , after passing the same events through a 
realistic GEANT detector simulation. The difference be- 
tween the two cases for the a, S and k are about an order 
of magnitude smaller than their absolute values. Typical 
values of such differences for central Au+Au 200 GeV 
collisions are -0.37±0.05, 0.02+0.05 and -0.06+0. 12 for 
a, S, and k, respectively. These results indicate that the 
effects on the shape of the distributions are small. The 
effect on the yields of p(p) is discussed elsewhere (il. Ilij. 
The systematic errors are estimated by varying the fol- 
lowing requirements for p(p) tracks: DCA, track quality 
reflected by the number of fit points used in track re- 
construction, and the dE/dx selection criteria for p(p) 
identification. The typical systematic errors are of the 
order 10% for M and a, 25% on S and 30% on k. The 
statistical and systematic (caps) errors are presented sep- 
arately in the figures. 

To understand the evolution of centrality dependence 
of moments in Fig. [21 we invoke the central limit the- 
orem (CLT) and consider the distribution at any given 
centrality i to be a superposition of several independent 
source distributions. We assume the average number of 
the sources for a given centrality to be equal to some 
number C times the corresponding (A part ), and obtain 
[21]: 

Mi - CM^V).. (1) 
af = Ca*(N palt ) h (2) 

ft = &/[^(Mil. ( 3 ) 

and 

Ki = Kx/[C{N pait )i\. (4) 

The various moments of the parent distribution M x , <t x , 
S x , k x and constant C have been determined from fits to 
data. The dashed lines in Fig. [5] show the expectations 
from the CLT. The x 2 /ndf between the CLT expecta- 
tions and data are < 1.5 for all the moments presented. 
If collision centrality reflects the system volume, then 
the results in Fig. [2] which approximate baryon number 
susceptibilities suggest that the susceptibilities do not 
change with the volume 0j. Deviations from (A par t) scal- 
ing could indicate new physics such as might result from 
the CP. 

To get a microscopic view, we present two observables, 
Sa and kct 2 , which can be used to search for the CP. 
These products will be constants as per the CLT and 
other likely non-CP scenarios, as seen from the depen- 
dences on (A/part) discussed above. These observables 
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are related to the ratio of baryon number susceptibili- 
ties (xb) at a given temperature (T) computed in QCD 



models as: Sa = 



X (2) /T2 



and 



Close to 



x w /7 » - - 

the CP, models predict the net-baryon number distribu- 
tions to be non-Gaussian and susceptibilities to diverge 
causing Sa and Ka 2 to deviate from being constants and 
have large values. Figure [3J shows that Sa and na 2 for 
Au+Au collisions at y/s^ = 19.6, 62.4, and 200 GeV 
are constants as a function of (A part ). 
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FIG. 3: (Color online) Centrality dependence of (a) Sa and 
(b) ko' 1 for AiVp in Au+Au collisions at \J snn = 19.6, 62.4, 
and 200 GeV compared to various model calculations. The 
shaded band for Sa and Ka 2 reflects contributions from the 
detector effects, (c) shows the model expectations for Ka 2 
from various physical effects in Au+Au collisions at 200 GeV. 
The lattice QCD results are for net-baryons corresponding to 
central collisions [91. See text for more details. 



In Fig.[3fa), lattice QCD results on Sa for net-baryons 
in central collisions are found to agree with the mea- 
surements. Near the CP, the system will deviate from 
equilibrium [l2| and results from lattice QCD, which as- 
sumes equilibrium, should not be consistent with the 
data. These lattice calculations, which predict a CP 
around fiB ~ 300 MeV, are carried out using two-flavor 
QCD with number of lattice sites in imaginary time to 
be 6 and mass of pion around 230 MeV [6j. The ratios 



of the non-linear susceptibilities at finite are obtained 
using Pade approximant resummations of the quark num- 
ber susceptibility series. The freeze-out parameters as a 
function of y / snn are taken from [22| and T c = 175 MeV. 

To understand the various non-CP physics background 
contribution to these observables, in Fig. [3] we also 
present the results for the net-proton distribution as 
a function of (A part ) from UrQMD (ver.2.3) [H], HI- 
JING Hi, AMPT (ver.1.11) 2|, and Therminator 
(ver.1.0) [25j models. The measurements are consistent 
with results from various non-CP models studied. In 
Fig. [UJc), several model calculations from Au+Au colli- 
sions at 200 GeV are presented to understand the effect 
of the following on our observable: with (W) and without 
(W/O) resonance decays, inclusion of all baryons (both 
studied using UrQMD), jet-production (HIJING), coales- 
cence mechanism of particle production (AMPT String 
Melting, ver.2.11), thermal particle production (Thermi- 
nator), rescattering (UrQMD and AMPT). All model cal- 
culations are done using default versions and with the 
same kinematic coverage as for data. The na 2 (Fig. [3b) 
and Sa (Fig. [HJi) are found to be constant for all the cases 
as a function of (A part ). This constant value can act as 
a baseline for the CP search. QCD model calculations 
with CP predict a non-monotonic dependence of these 
observables with (A part ) and v / snn (13j . 
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FIG. 4: (Color online) snn dependence of Ka 2 for net- 
proton distributions measured at RHIC. The results are 
compared to non-CP model calculations (slightly shifted in 
^/snn). The left-right arrow at the bottom indicates the en- 
ergy range for the CP search at RHIC. 

Figure @] shows the energy dependence of na 2 for AA p , 
compared to several model calculations that do not in- 
clude a CP. The experimental values plotted are average 
values for the centrality range studied; they are found 
to be consistent with unity. Also shown at the top of 
Fig. 2] are the fiB values corresponding to the various 
Jsnn 18, H3- We observe no non-monotonic depen- 
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dence with ^/snn- The results from non-CP models are 
constants as a function of ^/snn arid have values between 
1-2. The result from the thermal model is exactly unity. 
Within the ambit of the models studied, the observable 
changes little with change in non-CP physics (such as due 
to change in /ib, collective expansion and particle produc- 
tion) at the various energies studied. From comparisons 
to models and the lack of non-monotonic dependence of 
k<j 2 on j5nn studied, we conclude that there is no in- 
dication from our measurements for a CP in the region 
of the phase plane with < 200 MeV. It is difficult 
to rule out the existence of CP for the entire /^b region 
below 200 MeV. The extent to which these results can 
do that is guided by the following theoretical work. One 
QCD based model including a CP (£ = 3 fm) predicts the 
value of Kcr 2 to be at least a factor of 2 higher than the 
measurements presented (kct 2 ~ 2.5, 35, 3700 for the CP 
at v /ii5" = 200, 62.4, and 19.6 GeV, respectively) [H]. 
In addition, the expectation of the extent of the critical 
region in ^b is thought to be about 100 MeV frl l26j|. 

In summary, the first measurements of the higher 
moments of the net-proton distributions at midrapidity 
(| y |< 0.5) within 0.4 < p T < 0.8 GeV/c in Au+Au 
collisions at ^/snn = 19.6, 62.4, and 200 GeV have been 
presented. New observables Sa and na 2 derived from the 
AiVp distribution to search for the CP in heavy-ion col- 
lisions are discussed. These observables are found to be 
constant as a function of (iV part ) for all collisions energies 
studied. This is consistent with expectations from the 
central limit theorem and in general agreement with re- 
sults from various models without the CP. The measured 
Sa in central collisions are consistent with lattice QCD 
calculations of the ratio of third order to second order 
baryon number susceptibilities. Within the uncertain- 
ties, na 2 is found to be constant as a function of ^/snn 
studied. This trend is consistent with models without a 
CP and in sharp contrast to models [13] which include 
a CP in this /ib range. Our measurements show no evi- 
dence for a CP to be located at /iB values < 200 MeV in 
the QCD phase plane. The RHIC beam energy (100 < 
1-iB < 550 MeV) scan will look for non-monotonic vari- 
ation of na 2 for net-protons as a function of y / SNN to 
locate the CP. 
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